Three-photon energy-time entanglement 
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Entangled quantum particles have correlations stronger than those allowed by classical physics. 
These correlations are the focus of of the deepest issues in quantum mechanics 1-3 and are the 
basis of many quantum technologies. The entanglement of discrete particle properties has been 
studied extensively in the context of quantum computing ^, cryptography 5 , and quantum re- 
peaters [6] while entanglement between the continuous properties of particles may play a critical role 
in improving the sensitivity of gravitational wave detectors [T] |8], atomic clocks^, and other high 
^N^ precision instruments. The attributes of three or more entangled particles are fundamentally differ- 

^^ ent from those of two entangled particles [10f[14] . While the discrete variables of up to 14 ions[15j 

^^ and the continuous variables between three intense optical beams [IGJ [17] have been entangled, it 

fSI has remained an open challenge to entangle the continuous properties of more than two individual 

. particles. Here we experimentally demonstrate genuine tripartite continuous- variable entanglement 

^ between three separated particles. In our setup the three particles are photons created directly from 

^^ a single input photon; the creation process leads to quantum correlations between the colours, or 

^—i energies, and emission times of the photons. The entanglement between our three photons is the 

QQ three-party generalization of the Einstein- Podolsky- Rosen (EPR) 1 correlations for continuous vari- 

^Nl ables, and allows for new fundamental tests of quantum mechanics to be carried out. Our scheme 

can be extended to carry out multi-particle Franson interferometrv p^ 119] , and opens the possibil- 
f I ity of using additional degrees of freedom in our photons to simultaneously engineer discrete and 

1-5 continuous- variable hyper- ent angled states that could serve as a valuable resource in a wide variety 

of quantum information tasks. 
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Cu We directly generate three entangled photons using the nonlinear process of cascaded spontaneous parametric 

^ downconversion (C-SPDC)[20]. In downconversion, a pump photon, with frequency cjp, inside a nonlinear material 

, ^^ will occasionally fission into a pair of daughter photons with frequencies cjq and uji . The total energy in the process 
is conserved [2T] with huOp = hujQ + huji. The daughter photons share strong energy and time correlations that are 

^-H the hallmark of entanglement [HI [19]. The SPDC process is repeated with one of these daughter photons, at cjq, 

^ now serving as the pump, creating a pair of granddaughter photons simultaneously at uj2 and us- Again energy is 

^*^ conserved, and the total energy of the the three photons created in C-SPDC must sum to the energy of the pump: 

^p. hujp = hui + huj2 + hus. The simplified representation of our three-photon state in frequency space, assuming a 

vQ monochromatic pump, has the form 

O "^CSPDC ^ / / <icJidcj2Gi(a;i,a;p -a;i)G2(a;2,^p -^1 -^2)tti(^i)4(^2)4(^P -^1 -^2)|0), (1) 

where Gi{uoi^ujp — uJi) and 6^2(^2, ^p — uoi — UJ2) are the joint-spectral functions resulting from the phase-matching 
conditions of the first and second SPDC crystals respectively [22 . The three photons, consequently, share strong 
spectral correlations and exhibit genuine tripartite energy-time entanglement. 
^ To verify the tripartite entanglement of the photons generated in our C-SPDC process we use continuous variable 

Cd entanglement criteria, that we derive based on the work of van Loock and Furusawa[23]. for position and momentum. 
Consider three separable particles each described by the dimensionless observables x/^, Pk {k = 1,2,3) fulfilling the 
commutation relations [x^^Pi] = iSki (note: van Loock and Furusawa[23 use a different commutation relation than 
the one used here). Each individual particle must satisfy the uncertainty relationship AxiApi > 1/2. Together, all 
three particles must satisfy the following position-momentum uncertainty inequalities (see Supplementary Information 
for details): 

A(x2 - xi)A{pi +P2 +P3) > 1, (2) 

A(X3 -X2)A(pi +P2 +P3) > 1, (3) 

A(x3 - xi)A(pi +P2 +P3) > 1. (4) 

Violating any one of these inequalities is sufficient to demonstrate that a state contains some entanglement. Violating 
any two inequalities demonstrates that the state is fully inseparable [24 . For pure states full inseparability implies 
genuine tripartite entanglement [25^ . However, full inseparability and genuine tripartite entanglement are not, in 
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general, the same thing. Mixtures of bipartite entangled states that are fully inseparable but not genuinely tripartite 
entangled are also capable of violating two of the above inequalities. A more general entanglement criterion is therefore 
required to detect genuine tripartite entanglement. In the Supplementary Information we provide an overview of the 
definitions of full inseparability and genuine tripartite entanglement and derive the following inequalities: 

[A{x2 - xi) ^ A{x3 - xi)] A{pi ^ P2 ^ Ps) > 1. (5) 

[A(x2-Xi)+A(x3-X2)]A(pi+p2+P3) > 1, (6) 

[A{x3 - X2) ^ A{xs - xi)] A{pi ^ P2 ^ Ps) > 1, (7) 

[A{x2 - xi) ^ A{xs - xi) ^ A{xs - X2)] A{pi ^ P2 ^ Ps) > 2. (8) 

Violating any one of them is sufficient to demonstrate genuine tripartite entanglement. 

The position and momentum operators x and p are well-defined for narrow-band photons [2F, such as those generated 
by our C-SPDC process, with the usual commutation relation [x^p] = i. Because photons propagate at the speed of 
light, c, measuring the arrival time, t, of a photon at a single-photon detector is equivalent to measuring its longitudinal 
position X {t = x/c), and measuring its frequency, cj, is equivalent to measuring its longitudinal momentum p 
{huj = cp). Using this correspondence it is possible to write down the energy-time equivalents to the inequalities in 
equations |48l[52l 

[A{t2 - h) ^ A{ts - h)] A{uji ^ UJ2 ^ uJs) > 1. (9) 

[A{t2-ti)^A{t3-t2)]A{uJi^UJ2^UJ3) > 1, (10) 

[A{t3-t2)^A{t3-ti)]A{uJi^UJ2^UJ3) > 1, (11) 

[A{t2-h)^A{t3-h)^A{t3-t2)]A{uJi^uj2^uJ3) > 2. (12) 

States of the form in Equation [l] can violate all the inequalities maximally, ie the left-hand side goes to zero, and 
thus exhibit genuine tripartite entanglement. 

Measuring the difference in arrival times of the three photons using fast single-photon detectors gives the required 
timing uncertainties for testing the inequalities. However, directly measuring the frequencies of each individual photon 
with the precision needed (sub-GHz resolution over a bandwidth of several THz) to violate the inequalities is not 
technologically feasible. A solution to this problem is to use the fact that energy is conserved in the process of 
downconversion. The energy of the pump is equal to the energy of the three daughter photons created in C-SPDC 
{hwp = hwi + J1W2 + ^co's); measuring the frequency of the pump provides a direct measurement of the total frequency, 
of the three daughter photons required by the inequalities. This conservation of frequency has previously been tested 
using Franson interferometrvp^S } fT9 ] . and has been measured down to a line width of 200 kHz, an uncertainty much 
smaller than the scale considered in our experiment, in second-harmonic generation [^[28] (the time-reversed process 
of SPDC). 

To create our three entangled photons using C-SPDC (see Figure 1), first a narrowband pump laser at 404 nm 
is used to produce a pair of non-degenerate SPDC photons at 776 nm and 842 nm. The photon at 776 nm is then 
sent through a second SPDC crystal where a pair of granddaughter photons at 1530 nm and 1570 nm are generated 
(see Methods Summary for more details). This process leaves the 842 nm, 1530 nm and 1570 nm photons entangled 
in energy and time. Our setup detects an average of 7 triples/hour, from which we can infer the generation of 45 
triples/minute accounting for losses due to coupling and detection. In order to obtain sufficient photon counts with 
small statistical ffuctuations, data was collected for a total of 72.6 hours. 

The timing information from the detections was analyzed, and the triple coincidence counts binned into a 2D 
histogram based on ^2 — ti and ta — ^2 as shown in Figure 2. From the histogram it is clear that the photon 
arrivals are tightly correlated in time. The uncertainty in the arrival time between any pair of photons can be found 
by integrating over the arrival time of the other photon, removing its dependence as shown in Figure 3. From these 
integrated histograms we find that A{t2-ti) = 0.37±0.02 ns, A{t3-t2) = 0.162±0.004 ns, and A{t3-ti) = 0.31±0.02 
ns. Our measurements are limited by the timing jitter in our detectors and the resolution of our time-tagging unit 
(156 ps). The effect of the jitter can be clearly seen in the elliptical shape of the 2D arrival time histogram — the 
jitter on the detector used to detect the 842 nm photon is a factor of two larger than the jitter of the two telecom 
detectors. This is reffected in the uncertainty A(t3 — ^2) which is approximately a factor of two smaller than either 
A{t2 — h) or A{t3 — ti). Alternatively we can study the two-photon coincidences between detectors Dl & D2 and 
D2 & D3 independent of the third detector (see Supplementary Information) to verify that integrating over the third 
photon yields the correct two-photon timing histograms. The need for gating with the 1570 nm detector (^3) prevents 
the coincidences between ^3 and ti from being analyzed independent of ^2, but the 50 ns gate width is much larger 
than the uncertainty in the arrival time of a photon, and approximates the response of a free-running detector. 

The energy uncertainty of the photon triplets is given by the energy uncertainty in the 404 nm pump photons. To 
measure the uncertainty in the pump energy, a scanning Fabry-Perot interferometer (FPI) was used to continuously 
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FIG. 1: Our three entangled photons are created using C-SPDC. A narrowband pump laser at 404 nm downconverts into a pair 
of orthogonally polarised photons at 842 nm and 776 nm inside a periodically-poled Potassium Titanyl Phosphate (PPKTP) 
crystal. A filter (FP) removes the remaining pump light. A polarising beamsplitter is used to separate the two photons, and 
narrowband filters, FO and Fl, are used to block stray light. The photon at 842 nm is coupled into a single-mode fibre and 
sent to the single-photon detector Dl. The photon at 776 nm is coupled into single mode fibre and sent to a periodically-poled 
Lithium Niobate (PPLN) waveguide where it downconverts into a pair of photons at 1530 nm and 1570 nm. The photons 
are outcoupled into free space where a dichroic mirror is used to split the photons. The photons are then coupled back into 
single-mode fibre and sent to single-photon detectors D2 and D3 (see Methods Summary for more information about the 
detectors). The signals from all three detectors are sent to a time tagging unit, and a computer (PCI) is used to process 
coincidence events. The spectrum of the 404 nm pump laser is continuously monitored throughout the run using a Fabry- Perot 
interferometer (FPI) controlled by a second computer (PC2). 



monitor the bandwidth of the 404 nm laser throughout the experiment. Due to instabilities caused by temperature 
fluctuations, the measured bandwidth fluctuates over time as shown in Figure 4a, leading to the distribution in Figure 
4b. The average value and standard deviation of this distribution yield a pump bandwidth of Acc;p/27r = (6 ±2) MHz. 
The four measured time-bandwidth products for our three photons are 
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•A(t3-ti)]A(a;i+cj2+a;3) 
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0.03 ±0.01, 
0.02 ±0.01, 
0.018 ±0.005, 
0.03 ±0.01. 
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(16) 



Our three photons strongly violate inequalities [9p2] and are genuinely tripartite entangled. Our state exhibits 
energy-time correlations close to to the ideal state described in Equation [l] where the time-bandwidth products are 
exactly zero. This state is the continuous-variable analogue to the famous Greenberger, Home and Zeilinger (GHZ) 
entangled state[10, 23 , and the natural extension of the two-party continuous- variable EPR state [18, 19, 29 . 

Recent improvements in telecom wavelength detectors [30 and advances in nonlinear materials promise to dramat- 
ically increase our detected triples rate ^SJj i32] . Furthermore, new techniques to enhance the strength of nonlinear 
effects [33l [34] means that our scheme can in principle be scaled up to larger photon numbers. A major advantage 
of our states is that the continuous- variable entanglement is distributed amongst three individual photons, each at a 
different, tuneable, wavelength, enabling the creation of hyper-entangled states that simultaneously take advantage of 
both discrete and continuous variable quantum correlations. This multiplexing of entanglement over multiple discrete 
and continuous degrees of freedom may have important applications in quantum communication tasks. For example, 
a slight modification to our setup would enable a photon at 776 nm to be interfaced with an atomic storage medium 
like Kb while the remaining two photons are transmitted over telecom fibres to remote quantum nodes. This would 
open up new possibilities in the storage and distribution of quantum information needed for quantum computing, 
cryptography, and secret sharing, and could lead to new fundamental tests of quantum mechanics. 
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FIG. 2: 2D histogram of the timing information for the measured triple coincidences over 72.6 hours. The triple events are all 
localized to a small region of the histogram, indicating strong correlations in the arrival times of the three photons. 
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FIG. 3: Histograms of the difference in arrival times between two of the three photons measured over 72.6 hours. Each 
histogram is obtained by integrating the triples counts over the arrival time of the third photon to remove its dependence 
on the results. From this we find the uncertainty in the arrival times of the photons to be a) A(t2 — ti) = 0.37 ± 0.02 ns, 
b) A(t3 — t2) = 0.162 zb 0.004 ns, and c) A(t3 — ti) = 0.31 =b 0.02 ns. The timing uncertainties were verified using two-fold 
coincidence data that was obtained at the same time as the three-fold coincidence data (see Supplementary Information). 
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FIG. 4: Bandwidth of the pump photons as measured by a Fabry-Perot interferometer at five minute intervals during the 72.6 
hour run. a) The measured bandwidth as a function of time. The fluctuations in the measured bandwidths are the result of 
thermal drifts in the apparatus over the course of the run. b) Histogram of the measured pump bandwidths over the duration 
of the run. The average bandwidth is measured to be Ac(;p/27r =6 MHz with a standard deviation of 2 MHz (illustrated by the 
shaded regions in the graphs). 



A. Methods Summary 

In our setup (shown if Figure 1) we use a grating-stabihzed pump laser with a wavelength of 404 nm and a bandwidth 
of 5 MHz (Toptica Bluemode) to pump a 30 mm PPKTP crystal phase matched for Type-II SPDC. A pair of 
orthogonally polarised signal and idler photons at 842 nm and 776 nm respectively are generated co-linearly, and a 
polarising beam splitter (PBS) is used to separate them. The signal photon at 842 nm is coupled into an optical fibre 
and sent to a Si single-photon detector. With 12 mW of pump power 10^ signal photons/s are detected. The idler 
photon is fibre-coupled and sent to a second SPDC crystal, a 30 mm Type-I phase matched PPLN waveguide (HC 
Photonics), where it fissions into a pair of granddaughter photons at 1570 nm and 1530 nm. These granddaughter 
photons are outcoupled into free space and then split using a dichroic mirror. The photon at 1530 nm is sent to 
a free-running InGaAs/InP- Avalanche Photo Diodes (Princeton Lightwave, Negative Feedback Avalanche Diode - 
NFAD) detector cooled to 193 K operating at 10% efficiency with approximately 100 dark counts/s. This detector is 
used to gate a second InGaAs/InP detector (iD Quantique, id201-SMF-ULN) operating at 25% detection efficiency 
with a 50 ns gate window to detect the granddaughter photon at 1570 nm. The gated detector had a much higher 
dark count rate of approximately 5 x 10~^ dark counts/ (ns of gate). The arrival times of each of photons in the 
three detectors are recorded by a time-tagging system (DotFast/UQDevices) with 156 ps resolution. In this way all 
the timing statistics from the two-fold and three-fold coincidence events generated by the C-SPDC process can be 
measured. 
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I. SUPPLEMENTARY INFORMATION 

A. Fully inseparable versus genuine tripartite entanglement 

There are two definitions ttiat iiave been introduced for describing three-particle entanglement: fully inseparable 
[24] and genuine tripartite entanglement [25 . While these terms are sometimes used interchangeably in the literature, 
they mean different things. In [24 a classification scheme for three particle entanglement was developed based on the 
following set of tripartite states: 

P = ^ViPl,2,i^ P3,i^ (17) 

i 

P = X]^iPl,3,i(^p2,i, (18) 

i 

P = ^ViP2,34 ^ Pl,i, (19) 

i 

P = ^r]iPi,i^ P24^ P3,i- (20) 

i 

The first three states are biseparable as they are factorizable across a single cut while the final state is said to be fully 
separable as all three particles can be factorized. Tripartite states that cannot be written in any of these forms were 
defined to be fully inseparable [24] . 

Genuine tripartite entanglement [25 refers to states that cannot be written as a convex sum of just fully separable 
and biseparable states. In other words, they cannot be written as a convex sum of the states in equations [T7|[T8l 
These two definitions are not equivalent. For example, consider the state, 

P = 2^12 (^ P3 + ^Pl ^ P23' (21) 

This mixed state is not genuine tripartite entangled as it is a convex sum of biseparable states, yet it may be fully 
inseparable. In section D we present a more explicit example of a state that is fully inseparable but does not contain 
genuine tripartite entanglement. 

In [23 a criterion for detecting tripartite full inseparability in continuous systems is developed. In the following 
sections we build on this work to derive a stronger criteria capable of detecting genuine tripartite entanglement. 

B. Uncertainty relations 

Consider three separable particles each described by the dimensionless observables x^^ Pk (^ = 1^2, 3) fulfilling the 
commutation relations [xk^Pi] = iSki (note: the derivation in [23 uses a commutator with a different scaling factor). 
Each individual particle must satisfy the uncertainty relation AxiApi > 1/2. 

For a general three particle state p = ^- r]iPi23,i, the variance in A^(xi — X2)p is given as: 

A'(xi-X2)p = {{xi-X2f)^-{xi-X2)l (22) 



i \ i / i i 

i i \ i / 



(23) 

(24) 

(25) 
From the Cauchy-Schwarz inequality ^^ r]i {xi — ^2)^ > C^i ^i {xi — X2) ^Y ^ therefore the last two terms in equation 



25 will always be greater than or equal to zero. This leads to the inequality: 

A'(xi-X2)p > ^7/,A2(xi-X2)^. (26) 
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This result holds for any sum or difference of operators of this form. The uncertainty in {xi — X2) for the mixed state 
p is greater than or equal to the uncertainties in {xi — X2) of the components of the mixture. 

A^{xi-X3)p > Y,Vi^^i^i-^3)i, (27) 

A\x2-X3)p > ^riiA'^{x2-X3)i, (28) 

i 

A^{Pl+P2+P3)p > ^ViA^{Pl+P2+P3)i- (29) 

i 

These results lead to the following position and momentum inequality: 



A^Xi-X2)A^pi+p2+P3) > (^r?iA2(a;i-a;2)i W^r?J■A^(Pl+P2+P3)i 
Applying the Cauchy-Scharwz inequality to the right hand side of the previous equation yields: 

A'^{xi-X2)A'^{pi+P2+P3) > {^riiA{xi-X2)iA{pi+P2+P3)i\ , 



(30) 



(31) 



A{xi-X2)A{pi+p2+P3) > ^7]iA{xi - X2)iA{pi + P2 + P3)i. (32) 

i 

Following the same steps we also find: 

A(x2 -X3)A(pi +P2 +P3) > ^r/iA(x2 -X3)iA(pi +P2 +P3)i, (33) 

i 

A(xi -X3)A(pi +P2 +P3) > ^r]iA{xi-X3)iA{pi^p2^P3)i- (34) 

i 

C. Uncertainty relations for detecting continuous variable genuine tripartite entanglement 
In |23] van Loock and Furusawa study the permutations of a particular class of three-particle states defined as: 

P = ^ViPi,i2 (^Pi,3, (35) 

i 

where particles 1 and 2 can be entangled with one another, but are separable from particle 3, and derive a set of 
inequalities: 

[A{xi - X2)]^ + [A(pi +P2+ Ps)? > 2, (36) 

[A(a;2 - X3)f + [A(pi +p2+ Ps)]' > 2, (37) 

[A(a;i - X3)f + [A(pi +p2+ Ps)]' > 2. (38) 

It is important to note that in [23 a different scaling factor for the commutator is used which leads to the right-hand 
side of these inequalities having a value of one instead of two in that work. 

As the derivation of inequalities [36p8| depends only on the form of the commutator, they also hold for rescaled 



variables x^ = sx^^ p'j^ = Pk/s for any scaling factor s. Minimizing the expressions on the left of the inequalities with 
respect to the scaling factor s yields the following product inequalities: 

A(xi - X2)A{pi +P2 +P3) > 1, (39) 

A{x2-X3)A{pi +P2 +P3) > 1, (40) 

A(xi - xs)A{pi + P2 + P3) > 1. (41) 
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States that simultaneously violate any two of these inequalities are fully inseparable. Now consider the following p 
that is a mixture of all possible biseparable and fully separable states. 

P = ^ V^Pi^2 ^ Pi,3 + ^ Mift,13 ^ Pj,2 + ^ 7/cP/c,23 ^ Pk,l + ^ ^lpl,l ^ Pl,2 ^ Pl,3, (42) 

i j k I 

where rji^ jUj^ 7/c, and ui are probabilities with ^^T]i + ^j Pj + Xl/c 7fe + Xl^ ^^ ~ ^- With the correct choice of 
parameters it is possible for p to violate two of the van Loock and Furusawa inequalities simultaneously despite not 
containing any genuine tripartite entanglement (an explicit example is given in section D). A more stringent test is 
needed to look for genuine tripartite entanglement. 
Using the results from inequalities [32][34| we find that: 

A(xi -X2)A(pi +P2+P3) > ^r?iA(xi -X2)iA(pi +P2 +P3)i + ^MjA(xi -X2)jA(pi +p2+P3)j (43) 

i 3 

+ ^7/eA(xi -X2)feA(pi +P2 +P3)/c + ^^|A(xi - X2)lA{pi + ^2 +^3)/- 
k I 

If the component of the state pi^i2 is entangled across particles 1 and 2, then the lowest possible value of A(xi — 
X2)iA{pi +P2 +P3)i is zero. The remaining three terms are for cases where particles 1 and 2 are on opposite sides of 



a separable cut; from equation 39 the minimum value of the uncertainty products are A(xi — X2)jA{pi -\-p2 -\-P3)j = 
A(xi — X2)kA{pi -\- P2 -\- P3)k =~^{xi — X2)iA{pi -\- P2 -\- P3)i = 1- The uncertainty product in equation 43 can then 
be written as: 

A{xi-X2)A{pi^P2^P3) > ^^i(0)+ K^/i, +^7fc + ^^H (1), 

i \ J k I J 

j k I 

Following the same arguments it is possible to show that: 

A{xi-xs)A{pi^P2^P3) > ^Vi^^lk^^i^i- (45) 

i k I 

A{X2 - X3)A{pi + P2 + P3) > ^Vi + ^^ij + ^'^l■ (46) 



I 



Adding inequalities 44 and H5] yields the sum inequality: 



[A{xi - X2) + A{xi - xs)] A{pi + p2 + Ps) > l]^i + ^Mi + 2(^7fe + ^ 

* j \ k I 



VI ■ (47) 



The right-hand side is minimized by setting ^^ 7/^ = ^^ J^i = giving: 

[A{xi-X2)^A{xi-xs)]A{pi^P2^P3) > 1. (48) 

Similarly we find: 

[A{xi - X2) ^ A{x2 - xs)] A{pi ^ P2 ^ P3) > 1, (49) 

[A{x2 - xs) ^ A{xi - xs)] A{pi ^ P2 ^ P3) > 1. (50) 



Violating any of these sum inequalities [48pQl indicates that the state cannot be written with at most biseparable 



terms, and thus must be genuine tripartite entangled. 

It is also possible to add the three inequalities 44][46 together: 



[A{xi-X2)^A{xi-X3)^A{x2-X3)]A{pi^P2^P3) > 2 ^ T/i + ^ /i, + ^ 7fc + 3 I ^ ^H , (51) 

which is minimized when ^^ i^i = leading to: 

[A(xi - X2) + A(xi - xs) + A(x2 - xs)] A{pi + p2 + P3) > 2. (52) 

Violating this inequality also demonstrate the presence of genuine tripartite entanglement. 
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D. Example states 

Consider the following three states of three different particles written in the position basis: 

1 _(^j^Y -(^^^Y -(^^Y f ^1-^2 \^ 

'^1(^1,^2,^3) = -T^e v2-i; e V2-2; e V2-3; e v 2^, ; ^ (53) 

1 _r-LV -f^^Y -f^^Y / ^i-^3 \2 

^2(^1,^2,^3) = ^r^e V2-4; e V2-5; e v2-6; e v 2., j ^ (54) 

1 -T^^V -T^^V -T^^V / x2-.3 \2 

'^3(^1,^2,3:^3) = ^T^e V2-7; e V2-8; e V2-9; e v 2., j ^ (55) 

where A^i, A/'2, and A^3 are normalization constants and ai through aQ are adjustable width parameters. The width dc 
represents a correlation length for each of the states. When (Jc ^ two of the particles in each of the states become 
perfectly entangled with one another. The lower bound of the derived uncertainty relations 48](52 are saturated by 
pure states of the form shown in equations [53]|55l and are therefore tight. 

In general, full inseparability does not imply genuine tripartite entanglement. Consider the fully inseparable state 
composed of an equal mixture of ipi and ^^2 with ac = 0, (12,3,5,6 = 1, and cri,4 -^ 00. This mixed state, despite 
containing no component of genuine tripartite entanglement, simultaneously violate the van Loock and Furusawa 



inequalities [39] and [41] However, when tested against the new inequality derived in [48j it reaches a minimum value 
of a/2 > 1; our inequality is able to discern that this state contains no genuine tripartite entanglement. A similar 
fully inseparable mixed state, composed of an equal mixture of V^i, '^2, and i/js with ac = 0, 0-1,2,4,6,8,9 = 1, and 
^'"3,5,7 = 1/a/2, is capable of simultaneously violating all three of the van Loock and Furusawa inequalities despite 



not being genuinely tripartite entangled. When tested against the new inequality 52 it reaches a value of y/6 > 2, 
demonstrating that it is not genuinely tripartite entangled. 

In the specific case of pure states, full inseparability implies genuine tripartite entanglement. With the assumption 
of purity, both the van Loock and Furusawa inequalities as well as the ones derived here are capable of detecting 
genuine tripartite entanglement. For example, the pure state 

1 _(^J^Y -(^^Y -(^^^Y / a3i-cc2 \2 / a;i-cc3 \2 

tlj4{xi,X2,xs) = -—e v2-i; e V2-2; e V2-3; e v 2.^ J e v 2ae J ^ (56) 

where N4 is a normalization constant, exhibits perfect correlations between particles 1-3 when (Jc ^ (for finite 
widths cTi 2,3), and is both fully inseparable and genuinely tripartite entangled. This pure state will violate the van 
Loock and Furusawa inequalities as well as the inequalities [48pQ] and |52J It is only in the more general case of mixed 



states that the relations derived by van Loock and Furusawa need to be extended in order to detect genuine tripartite 
entanglement. 

II. ALTERNATE MEASUREMENT OF THE PHOTON TIMING UNCERTAINTY 

To verify that integrating over the third photon yields the correct two-photon timing histograms, the raw two- 
photon coincidence data from the measured time tags between photons 1 and 2 as well as the coincidences between 
photons 2 and 3 were analyzed independently of detecting the other photon. While the timing uncertainty between 
two detectors in Equations 2-4 are conditioned on the presence of a third photon, this measurement of the arrival 
times (shown in Figure [5| is different in that it is conditioned only on two photons. From this data, it was found that 
A(t2 — ti) = 0.4±0.2 ns and A(t3 — ^2) = 0.16±0.04 ns, which agree with the integrated values measured in the paper. 
The reason for the larger error is due to the large number of accidental counts. In the histogram between ^2 — ti the 
large background is caused by the coincidences between dark counts on the free-running InGaAs/InP detector and 
the 10^ trigger photons per second on the Si detector. The number of dark counts, and hence signal-to-noise ratio, is 
lower in the ^3 — ^2 histogram; here the background is primarily due to accidental coincidences between dark counts 
between the two telecom detectors. As the first telecom detector ^2 gated the second telecom detector ts, our setup 
did not allow us to measure the two- fold coincidence histogram ts — ti independent of photon 2. 
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FIG. 5: Histograms of the arrival times between a) photons 1 and 2 and b) photons 2 and 3. This was accompUshed by examining 
the two-fold coincidences, as opposed to the three-fold coincidences, from the collected time tags. From these histograms an 
uncertainty in the difference between the arrival times was found to be A(t2 — ti) = 0.4 zb 0.2 ns and A(t3 — t2) = 0.16 =b 0.04 
ns which agrees with the integrated uncertainties measured in the main paper. 
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FIG. 6: Direct measurement of two-photon energy-time entanglement. A narrowband pump laser at 404 nm is used to pump 
a periodically-poled Potassium Titanyl Phosphate (PPKTP) that downconverts into a pair of orthogonally polarized photons 
at 842 nm and 776 nm. A filter (FP) removes the remaining pump. A polarizing beamsplitter is used to separate the two 
photons, and narrowband filters, FO and Fl, are used to minimize stray light. The photons at 842 nm and 776 nm are coupled 
into single- mode fibres and sent to the single-photon detectors DO and Dl. The signals from both detectors are sent to a 
time tagging unit, and a computer (PCI) is used to process coincidence events. The spectrum of the 404 nm pump laser is 
continuously monitored throughout the run using a Fabry-Perot interferometer (FPI) controlled by a second computer (PC2). 



III. TWO PHOTON ENERGY-TIME ENTANGLEMENT 



In a separate experiment we directly measured for 1.6 s the two photon energy-time entanglement between the pair 
of daughter photons (at 842 nm and 776 nm) created at the first stage of our C-SPDC process as shown in Figur^ 
Both daughter photons were sent to single-photon detectors. The pump power was attenuated so that approximately 
14,000 coincidences/s were detected. A time tagging unit was used to store the arrival times of the photons, and 
a computer then sorted the time tags into extract coincidence events. The histogram of arrival times of the two 
photons is shown in Figure [7| and the uncertainty in the arrival times measured to be A{to — ti) = 0.30 ± 0.01 
ns. Simultaneously the Fabry-Perot interferometer monitored the bandwidth of the pump, yielding a bandwidth of 
Aujp = 4.6 ± 0.8 MHz. 

To show bipartite entanglement it is sufficient to violate the two party inequality A(co'o -\- oJi)A{to — ti) > 1 [23] , 
As energy is conserved in SPDC we can use the fact that ujq -\- Ui = ujp. This leads to a time-bandwidth product of 
A(a;o + uJi)A{to — ti) = 0.0014 ± 0.0002 that violates the classical limit by more than 4000 standard deviations. 
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FIG. 7: Two-photon coincidence histogram of the difference in arrival times of the 842 nm and 776 nm photons from downcon- 
version. From this histogram, the uncertainty in the difference in arrival times was measured to be A (to — ti) = 0.30 it 0.01 
ns 



